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APPLICATION OF THE METHOD OF STEEPEST DESCENT TO LAMINATED SHIELD 
WEIGHT OPTIMIZATION WITH SEVERAL CONSTRAINTS - THEORY 

by Gerald P. Lahti 
Lewis Research Center 

SUMMARY 

The method of steepest descent used in optimizing one -dimensional layered radiation 
shields has been extended to multidimensional, multiconstraint situations. The multi- 
dimensional optimization algorithm and equations have been developed for the case of a 
dose constraint in any one direction being dependent only on the shield thicknesses in that 
direction and independent of shield thicknesses in other directions. Expressions are de- 
rived for one-, two-, and three-dimensional cases (one, two, and three constraints). 

The procedure is applicable to the optimization of shields where there are different dose 
constraints and layering arrangements in the principal directions. 


INTRODUCTION 

One -dimensional shield weight optimization procedures are currently in use for the 
design of layered radiation shields (refs. 1 to 4). In general, though, shields such as for 
space power reactor applications may have different dose constraints in the several dif- 
ferent directions. This requires either the superposition of two (or more) separate one- 
dimensional optimizations or development of a procedure which optimizes in two or more 
directions simultaneously. The basis for a computer code to do the latter is developed 
in this report. Expressions for the general case of constraints being functions of all 
shield layer thicknesses are derived, but the complete algorithm for the optimization 
procedure is made in this report only for the case where the dose constraint in each 
principal direction is a function of thicknesses in that direction only. This has applica- 
tion to the optimization of shields where there are different dose constraints and layering 
arrangements in the different principal directions and there is a weak (or no) dependence 
of the dose in one direction on the thickness in another direction. 



The development of the method described in this report is parallel to the one- 
dimensional method developed by Bernick (ref. 3), programmed as the OP EX code 
(ref. 4), and later reprogrammed and reinterpreted as the OPEX-II code (ref. 5). 


MATHEMATICAL STATEMENT OF THE PROBLEM 

Assume a shield consisting of n layers. (For purposes of mathematical develop- 
ment, only the total number of layers need be considered. The n layers maybe con- 
sidered to be divided into ^ layers in the first direction, layers in the second di- 
rection, and so forth, with n = ^ + n 2 + . . . . ) is the thickness of the ith layer. 

Let a vector x be defined in n-dimensional vector space with quantities x^, Xg, 

. . . , x n as components : 



Let the total shield weight be denoted by W(x). The function W(x) will in general be 
nonlinear . 

The problem is then to minimize total shield weight subject to m dose constraints, 
say in certain directions. Mathematically, this can be expressed in the following way: 
Let 


W(x) — minimum 


with constraints 


D° = Dj(x) - C x = 0 


(1) 


(2a) 


2 



(2b) 


D° = D 2 (x) - C 2 = 0 



= D (x) - C = 0 
m ' m 


(2m) 


and 


x. a 0 (3) 

ii. 

where Dj(x) is defined as the value of the dose rate at the j point (i. e. , in the j 

direction) and is, in general, some nonlinear function of x, and is the required dose 

rate at point j. (Again, the x. and hence D. may be considered to be subdivided into 

1 J fh 

the different directions with D. strongly dependent on those x^ in the j direction and 
weakly dependent on, or completely independent of, the thicknesses in the other direc- 
tions. The present derivation first considers the general case of Dj dependent on all 
Xj.) 

Constraint (3) is necessary to ensure a solution with all Xj & 0, that is, a feasible 
solution. 


THE ALGORITHM 

The method of steepest descent as applied in this report is an iterative procedure to 
find a vector x m - n which is a solution to the problem stated in equations (1), (2), and 
(3). The procedure is as follows: 

(1) Choose an initial feasible vector Xq, (i. e. , one which satisfies constraints (2) 
and (3)). 

(2) Find a unit vector u, at point Xq, which points in the direction of maximum al- 
lowable decrease in the function W(x). The direction of u is chosen subject to the con- 
straint that u is tangent to the constraint surface defined by D°(x) =0 for all j (i. e. , 
the dose constraint is met). 

(3) Calculate a vector 


x l = 


Xq + X 0 U 


where X° 


is some interation parameter, a constant, say. 


3 


(4) Calculate D°(x^) for all j. Because of the nonlinear nature of D, D?(xj) will 
not, in general, be equal to zero. A corrected point x^ is then determined such that 

D?^) = 0 

(5) Set x x — Xq; repeat steps (2) to (5) until some convergence criterion on weight is 
satisfied. The final value of x thus computed should then be a close approximation of 

x min' 

DERIVATION OF u 
General Considerations 

Let vectors g, and a j, a 2 , . . . , a m be defined by 

g = VW(x) 
a x - vrfj(x) 

a 2 = VD°(x) 

a = VD° (x) 
m 

where 

' _ 8 _ 

0Xj 

8 

8x 2 

V = 

8 

. 9x n 



4 



and V is the conventional gradient operator as applied in n-dimensional space. The 
problem is then to find a vector u at some point x such that 


g . u — minimum 

(4) 

u • u = 1 

(5) 

aj • il = 0 

(6a) 

a 2 • u = 0 

(6b) 

a: • u = 0 

m 

(6m) 


Expression (4) indicates that the component of u along the direction of greatest increase 
of W is to be minimized (i. e. , the component in the direction of greatest decrease of 
W is to be maximized). In the case of no dose constraints, for example, u and g 
would point in opposite directions or 


g • u _ 
g 

Expression (5) constrains u to be a unit vector. This is added for computational con- 
venience. Expressions (6) simply state that u must be perpendicular to VD? (i. e. , u 
must lie in the plane of, or in general only be tangent to, a hypersurface of constant 
dose, D? = 0, for all j). 

Equations (4) to (6) can be solved by the method of Lagrange multipliers. Write the 
Lagrangian as 


!£= g ■ u + • u) + a 2 (a g • u) + . . . -j3((u • u) - 1) (7) 

where a j, a 2 , . . . Qf m , and 0 are undetermined multipliers. The stationary points 
of the objective function g • u are desired at a particular point x. At this particular 
point, g is fixed and only the components of U are variable; that is, g • u = g • u(uj) 
only. Therefore, to solve for the unknown u^ and multipliers on and /3, each partial 
derivative of equation (7), with respect to u^ is set equal to zero to determine sta- 
tionary points. That is, 


5 



^=0 

d<£ 
du „ 


= 0 


ay 


= 0 


or simply 




Evaluating this term by term, where Uj Is a unit vector in the i th direction, 


n 



n 


Yj s A =s 


i=l 


Similarly, 


V“i a i 


u) = 


“i a i 


and 


6 



V s ' 



U 1 + "2 ' u 2 + • 


+ u 

n 




= 2u 


or finally 


V_if = 0 = g + q-jSj + a 2 a 2 + . . . - 2/3u (8) 

Derivation of Oj 

Equation (8) along with constraints (5) and (6a) to (6m) provide a set of m + 2 equa- 
tions in m + 2 unknowns (u, /3, a^, ar 2 , . . . , ar m ). To solve this set of equations, 

take a^ • [eq. (8) ] and obtain 

g • aj + a-j a^ + o^aj • a 2 + . . . + o^aj • a m - 2/3aj • u = 0 (9a) 

Then take a 2 ■ [eq. (8) ] and obtain 

g • a 2 + a^ 2 • a 1 + ctf\ + . . . + • a m - 2/3a 2 • u = 0 (9b) 


m 


at* a 
1 m 


a l + 


<*2 a m 


a 2 + 


a a 
mm 


m 


2/3a 


m 


u = 0 


(9m) 


Because of boundary condition (6), terms like a^ • u vanish. Hence, equations (9) com- 
pose a set of m equations in m unknowns, namely, 


7 



( 10 ) 


' a l 

a l a l ' 

a 2 ‘ • 

a l* 

3-™ 

m 

'“l ‘ 


g • 

a l] 

a 2 

a l a 2 ‘ 

a 2 * - * 

a 2 * 

a^ 

m 

a 2 

= 

-S’ 

a 2 

_ a m ‘ 

a l a m ’ 

a 2 • • • 

a m’ 

a. 

“m. 


- g • 

3 

m 


In principle, equation (10) can be solved for any m. Consider the following cases 
m = 1: 




( 11 ) 


m = 2: 
Let 


D 


2 



a-i a i • a- 2 \ 

aj a 2 • a 2 J = ( *1 ' a l )(a 2 ' a 2> " (a 2 * a l><*l * a 2> 


Then 


a l a l ' a 2 


“1 = 


-g • a 2 a 2 • a 2 -(g • • a 2 ) + (g • a^ • a 2 ) 


D c 


(aj • a x )(a 2 • a g ) - (a g - a^ 2 


(12a) 


ra^ • a^ -g • a^> 




f 2 • aj -g • a 2 / -(aj • a x )(g • a g ) + (g • a 1 )(a 2 • a x ) 


D 0 


(aj • aj)^ 2 • a 2 ) - (a 2 • a^ 2 


(12b) 


8 



m = 3: 


Let 



a 2 a l 
a 2 a 2 



( a l • a ^)(&2 ' a 2)( a 3 ‘ a 3 ) 

-( a l • a j)( a 2 ' a 3) ~ ( a j ‘ a 2 ) ( a 3 

+(a j • a 2 )(a 2 ■ a^aj • a 3 ) 

+(a^ ■ ag) (a 2 • a j)(a 3 ■ a 2 ) 

-(a t • a 3 ) 2 (a 2 • a 2 ) 


a 3> 


Then 


a 


1 


/-g ‘ a. 


-g • a. 


t-g 




(13a) 


a 


2 



-g ’ a x 

-g • a 2 

-g • a 3 



°3 


a 


3 



-g • aA 

-g • a 2 j 

-g • y 


(13b) 


(13c) 


9 



Derivation of p and u 


To solve for /3, rearrange equation (8): 


u = 


g + + a 2 a 2 + * ' 


a a 
• m m 


2/3 


Take u- of both sides. Using both equation (14) and u • u = 1 results in 


u • u = 1 


(g + a.a, + a 2 a„ + . . . + a m ) (g + a 1 a 1 + « 2 a 2 + . . . « m a m ) 


2/3 


2/3 


One constraint (m = 1) . - For m = 1, equation (15) reduces to 

4/3 2 - g • g + 20-^ • g + a^aj • a^ 

From equation (11) 


O' 


-g • a J 
1 5 1 • *1 


So 


4/3 2 = g ■ g - 2(3 — -M (a x • g) +(- ^ ) (a x • = g * g - — 


f g ■ a 


1 \ /*- 




(g • a x ) 2 


V a l • a i/ 


a i ‘ a iy 


a i* a i 


4- ~ 

± yg • ~ 


2/3 = ± vg • g - 


a i ' a i 


The sign of this has not yet been resolved. 

From equation (14) and the value of a ^ from equation (11) we have for m = : 


( 14 ) 


(15) 


, then, 


10 



u = 


(16) 



For this to be a real solution, it is necessary that 


(g • a^ 21 
g • g 


a l ' a lJ 


>0 


(17) 


Finally, the sign must be resolved. 

Recall the original problem was to find the minimum u • g. Therefore, taking 
g- [eq. (16)] yields 


u • g = 



• g) 



g - 



(18) 


Or, with 


(g • a i) 2 
B = g • g - — — 

a l- a l 


Equation (18) becomes 


u • g = 


B 


±B 1//2 


Because it is necessary that B > 0 for a solution, the plus sign maximizes u • g and 
the minus sign minimizes u • g. Therefore, the minus sign is the required solution for 
this problem. 


11 


Hence, the final solution for the m = 1 case is 


g • a 


-g + 


u = 


ia l ' a l; 


1 V 


21 l/ 2 

(g ■ a,) 27 

g • g 


( 19 ) 


(a-. • a,) 


l'J 


This is the one-dimensional result derived by Bernick (ref. 3). 

Two or three constraints (m = 2 or m = 3). - From equation (15), for m = 2, 

4/3 = (g + + a> 2 % 2 ) ■ (g + + a^) 

„ „ _ 2_ _ _ 2_ _ 

= g • g + 2a^g • aj + 2a 2 g • a 2 + a-ja^. • aj + 2of 1 Q' 2 a 1 • a g + « 2 a 2 • a 2 (20) 


This may be expanded by using the appropriate expressions for a, equation (12). From 
equation (14), 


g + o' 1 a 1 + Q; 2 a 2 

u = 2/3 (21) 

where a ^ and q? 2 are obtained from equation (12) and /3 is obtained from equation (20). 
This is the general expression for the vector u for a two-constraint problem. 

In general, the dose in any direction is known and may be expressed approximately 
as a function of thicknesses in that direction only. This report presently generates the 
optimization algorithm for this special case only. That is, consider the dose at points 
in the first and second (and third) directions determined only by the thicknesses in that 
particular direction. It is stated previously in this report that the n thicknesses gen- 
erally might be regarded as being nj thicknesses in the first direction, n 2 in the sec- 
ond, and n^ in the third and that ^ + n 2 + n g = n. If the dose equations are 

D° = D°(xj) i = 1, nj 

D 2 = D 2 (xj) i = (nj + 1), (nj + n 2 ) 

Dg = Dg(xj) i = (n x + n 2 + 1), n 


12 



one obtains 


3D° 
3x j 


0 


0 

3D° 

3x 2 


0 


0 

3D° 


0 


o 

3x 

n l 



u 

0 


3D 2 


0 


3x„ 

n^+1 


0 


3D° 


0 


*2 = 

8x n 1 H-2 

a 3 ~ 

0 


0D° 


0 


Sx ni + n 2 


0 


0 


9D S 




9x n 1 +n 2 +l 

0 


0 


8 D“ 




0X n in n 

n j+ng+2 

0 


0 


3D^ 

9x n 


13 










As a result of this dose model, all terms of the form a^ • a. equal 0 for i / j. Al- 
though the dose terms in the different directions are decoupled, the weight derivatives 
are not. 

For the two-constraint, two-direction problem (m = 2), the values of a, equa- 
tion (12), reduce to 


-g' a l 

• a 


and 


a 


2 


-g ‘ a 2 
a 2 ' a 2 


From equation (15) 


4/3 


2 ^ ^ 


g • g + 2of 1 a 1 • g + 2a 2 a 2 


„2„2 


.22 


g + tfja^ + &2 a 2 


(a.. • g) 2 (a ? • g) 2 (g • a.) 2 (g • a ) 2 

= g- g - 2 _ 1 - _ - 2-4 — - + _ 1-+— ~ 

a l a l a 2 a 2 a l a l a 2 a 2 


= g • g - 


(aj ' g)‘ 


a r a i 


■ a 2 
a 2 a 2 


20 



g - 


(a- 1 ■ g) 2 
a l ' a l 


^2 

a 2 



Finally, from equation (14) and the expressions for 


Qfj, Q! 2 , and 0 just given, 


14 



( 22 ) 


g * a, g • a„ 

5 i 


U = 


a l ' a l 


a 2 ’ a 2 


(a< * g) (a» • D' 

g • g - 


, 1/2 


a l ' a l a 2 a 2 J 


For the reasons given previously, the minus sign in the denominator is again selected. 
Similarly, for m = 3, from equation (13), 

E >3 = (a t • ajXag • a 2 )(a 3 ■ a g ) 


a l = 


(-g • a 1 )(a 2 • a 2 )(a 3 • ffg) -g • a.^ 


Dr 


a i’ ^ 


a 2 ~ 


(-g • a 2 )(aj • aj)(a 3 • a g ) -g • a^ 




a 2 * a 2 


a 3 = 


(-g • a 3 )(a 2 • ag)^ • a^ -g ■ a. 


Do 


a 3 ’ a 3 


From equation (15), 



= g ‘ 


g • 


o _ 2 2— — — — — — — — 

g + aj( a 1 ■ a x ) + a 2 (a 2 • a 2 ) + a 3 (a 3 • a g ) + 2a^{g • aj) + 2a- 2 (g • a g ) + 2Q- 2 (g • a g ) 

(g • aj) 2 (g ‘ a 2 ) 2 (g-a 3 ) 2 2(gL- a^ 2 2(g • a 2 ) 2 2(g • a g ) 2 

(a^ • a^) (a 2 • a g ) (a g • a g ) (a^ • a^) (a 2 • a 2 ) (a g ■ a g ) 

(g • a x ) 2 (g ■ a 2 ) 2 (g • a 3 ) 2 

( a ! • a l> < a 2 • a 2 } (a 3 ‘ a 3 } 


15 



From equation (14), 


g • a 1 g • a„ 

g " I— ^ a l _ = a 2 


g 


u = 


a l ‘ a l 


a 2 ‘ a 2 


a 


~ ^3 


3 ' 3 


(g • a x ) 2 (g • a 2 ) 2 (g • ag)* ^ 

g . g 


a l ‘ a l 


a 2 ' a 2 


a 3 ' a 3j 


(23) 


Again the minus sign in the denominator is selected. 

By comparing equations (19), (22), and (23), one could extend the derivation of u to 
even higher dimensions (m > 3). This is, however, only of academic interest because 
one could probably not define a dose -thickness relation to fit the m-dimensional case. 


COMPLETION OF THE OPTIMIZATION ALGORITHM 


Let the vector u, evaluated by equation (19), (22), or (23) at x, , where k indicates 
the k tn optimization step, be denoted by u^. After u k and are determined from 

x k+i = x k + X\ (24) 


a dose rate calculation is made. If the quantity D°(x^ + ^ 
ferent from zero, a correction to x^+i must be made, 
ponents dx^, . . . , dx n is sought such that 


) (see eq. (2)) is sufficiently dif- 
An increment dx k+ j with com- 


- 5 «% + i> = 0 


(25) 


where 


®<W = VD X + i> • dx k + i 


(26) 


Also a minimum |dx k+ ^ | is desired. 

This problem is then to minimize dx ■ dx with the constraint 


6D(x k+ i) = D“(x k+ ,) 


16 



Again define a Lagrangian if by 


y = dx • dx + r[6D(x k+1 ) - D?(x k+1 )] 


(27) 


where y is an undetermined multiplier. We wish to find stationary points of the objec- 
tive function dx • dx at x k+1 . To solve for the unknown dxj and multiplier y, at the 
stationary points, again take derivatives of if, equation (27), with respect to dx. and 
set each equal to zero. The result is 


= 0 = 2dx 1 + yVDf(* +1 ) ■ dxj 




Sdx 


— = 0 = 2dx + yVD?(x k+1 ) - dx 
adx„ n 1 lc+i 71 j 


(28) 


n 


where dxj is a unit vector in the direction dx^. Equations (28) may be written as 


0 = 2dx + 7 VD°(x k+1 ) 


(29) 


Equations (25) and (29) provide a set of n + 1 equations in n + 1 unknowns dx^ and 
y. To solve this, from equation (29), we have 


dx = -X VD?(x k+1 ) 


Substituting this in equation (25) and using equation (26) give 


D ?<*k + 1> - VD°(x k + 1 ) ■ (2) VD?(x k+1 ) - 0 


or 


- 2D ?^k + l) 

v 2o °(x k+ i) 


17 



It follows, then, that 


^ _ D i (x k+1> VD f ( *k+l) 
V 2n “(x k+1 ) 


nO,- 


,th 


The new corrected vector x^^ which satisfies the i dose constraint is given by 


x k+l _ x k+l 


+ dx 


The procedure is to be repeated for each of the j directions. 

The algorithm is terminated when the relative weight change from one iteration to 
the next is less than some prescribed e. That is, when 


W(\+l> - W(x k ) 

w<*y 


< £ 


CONCLUDING REMARKS 

The algorithm for optimizing a two- or three-dimensional shield has been developed 
by extending a one -dimensional algorithm. Equations have been derived which permit 
implementation of this algorithm in a computer program similar to existing one- 
dimensional optimization codes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 17, 1971, 

112-27. 
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